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^ ■ Abstract 

In this paper, we investigate the scattering of two t/> 6 kinks and derive the 

> 

I/-) | real dynamics by solving the appropriate field equation numerically employing 

IT) '. 

. a Runga-Kutta method. We also use a collective coordinate approximation to 

find approximate dynamics, with the objective being to compare the approx- 
CO . imate motion to the real dynamics in order to test the validity of this kind of 



approximation, which is used extensively in the study of solitons. 



^- 1 Introduction 

Topological solitons are smooth, stable, particle-like solutions to non-linear wave 
equations that occur in field theories [I]. These include solitons such as monopoles in 
(3+1) dimensions and vortices and lumps in (2+1) dimensions. However, although 
some of the physically interesting theories with solitons may be exactly solvable 
in the static case (through the reduction of second order to first order differential 
equations via a Bogomolny argument [2]), they are usually not integrable upon the 
introduction of time dependence. This makes it difficult to study soliton dynamics 
and usually one must make use of approximations or numerical simulations. Here, 
we present both of these approaches, this time in the context of 6 kinks. Our aim 
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is to use one method to test the other. Section 2 presents a discussion of <fi 6 kinks. 
Then section 3 introduces our collective coordinate approximation, with section 3.1 
discussing one 6 kink and section 3.2 focusing on two kinks. Section 4 presents 
dynamics for two kinks. The paper ends with a conclusion. 

' kinks 

The 6 model is a field theory in (1+1) dimensions involving a real field (p(t, x). The 
Lagrangian density is given by 



where 



U{cf ) )= l - ( t>\l-ff. (2) 



Varying (TjQ) leads to the equation of motion 

d 2 6 d 2 (b dU 



dt 2 dx 2 

The Bogomolny equations are 



+ TJT = 0- (3) 



/>' = ±V2U. (4) 



For the 6 model they become 



0' = ±0(l-0 2 ). (5) 

This ODE can be solved exactly for both signs. Taking the negative sign, we find 
a kink (with 0(— oo) = — 1 and (f>(oo) = 0) and an antikink (with <f>(—oo) = and 
4>{oo) = —1). The kink solution is 

4>\(x) = . 1 , (6) 

y/l + exp (2(x + a) + log3) 

where a has been introduced as a constant of integration. Note that the factor log 3 

has been included to position the kink at the origin when a = 0. 



1 Note that kink- antikink scattering in the 4> e model has also been investigated recently, see [3] 
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Now taking the positive sign in (jSJ), we again find a kink (this time with <p(— oo) = 
and </>(oo) = 1) and an antikink (with <p(— oo) = 1 and <j)(oo) = 0). The kink 
solution is 



<h(*) = 1 =, (7) 

v/l + exp(-2(x-a)+log3) V ; 

where a is again an arbitrary constant and log 3 is included. 

In both ((6]) and ((7j), the constant a represents the position or modulus of the 

soliton. <j>i and <p2 represent kinks positioned at —a and a, respectively. Note we 

have the symmetry 

M-x) = -Mx)- (8) 

We can demonstrate that —a and a are the moduli of the kinks by considering the 
energy density 

£= l -^ + U = ^. (9) 
Here is a static kink field. (Q implies 

f- = Wtf. (10) 
ax 

Consider ()6]). When x = —a, <p\ ^ 0. Hence, for this point to represent the max- 
imum of the energy density and therefore the kink position we need <f)"(—a) = 0. 
Differentiating ([6]) we see that this is indeed the case. Similarly we find 4>2 ( a ) = 0- 
The dynamics of a kink given by ([6]) or (J7J) can be found by a simple Lorentz 
boost, e.g. 

fa(t,x) = -= ~ l (11) 

'l + exp (2{ 1 (x-vt + a+^))) 



where v is, as usual, the kink speed (with — 1 < v < 1) and 7 = ^ ^ is the 
Lorentz factor. Note that while the dynamics of one kink, in this model, can be 
found by direct integration, those of multi- kinks cannot. This is a result of the fact 
that the (J) 6 model is not an example of an integrable system and therefore, unlike 
for the well-studied sine-Gordon model, solution generating techniques such as the 
Backlund transformation are not applicable. 



3 



3 The collective coordinate approximation 



3.1 One kink 

In [1], a moduli space approximation was introduced. This kind of approximation 
was originally implemented for a system of solitons for which there are no static 
forces. In the case of monopoles it is possible to have multi-soliton configurations 
that are static with no forces between solitons (critical coupling) because of a precise 
balancing between a magnetic force and a scalar interaction (due to a Higgs field). 
This leads to geodesic trajectories of multi-solitons on the space of static configura- 
tions (the moduli space). For kinks, of course, only a single kink can be static and 
the calculation is simple. We demonstrate this for a 6 kink. 

A static kink is given by (J7J). Here the parameter or modulus is a real number 
a and hence the moduli space is simply M. The moduli space approximation is 
to consider a flow on this manifold, i.e. to make a time-dependent, and to then 
substitute this dynamical field into (0Q). Integrating over all space gives the dynamics 

L=-a 2 --, (12) 
8 4' v ; 

hence motion at a constant speed. Generally, when there are no static forces, there 

is geodesic motion on the moduli space. 



3.2 Two kinks 

When there are static forces, the moduli space approximation becomes harder to 
implement. We consider multi-solitons that experience static forces. Here, it is 
possible to make progress by taking 1-soliton configurations that are static and 
patching these into one field then performing a flow. When the solitons are far 
apart the field will be accurate. It is then hoped that when the solitons approach 
and interact, the approximation will still be good. This approach has been applied 
for solitonic structures such as lumps, vortices and monopoles, as well as kinks. See, 
for example, [HI [6] for a study of classical and quantum monopole scattering, [7J [8] 
for vortex scattering, and [9] for classical lump scattering. There is confidence in 
the validity of such a method. To provide further confidence we can perform "tests" 
on one-dimensional models. Such a test was performed on the scattering of two 
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sine-Gordon kinks by Sutcliffe [TO]. Here we show how to perform such a test on 
the (p 6 model. We choose to obtain a 2-kink field by patching together two 1-kink 
solutions in the following way 

(13) 



v /l + exp(2(x + a)+log3) ^/T+ exp (— 2(x — a) + log 3) 
where a G M + is the collective coordinate, so that Ai is one dimensional. To perform 
the truncation to a one dimensional system we substitute the ansatz (1131) into the 
Lagrangian density (TT1), and integrate over all x. This leads to the Lagrangian 



L = ^g(a)a 2 -V(a). (14) 



Here g(a) is the metric on Ai 



. . 1 f°° e ia dx , . 

g( a ) = - - 18 / 3 ~- 15 

2 ( 1+ 3 e -2(x-a))|( 1 + 3 e 2(x+a))| 

and V(a) is the potential 

/oo 
(9Vf + V 2 2 (l - V 2 2 ) 2 )dx, (16) 
-oo 

where 

( p 2{x+a) -2(x-a) 
~ 3+ ~ 5-1 (17) 
(l + 3e 2( a; +a))| (i + 3e- 2 (- a ))f ' 

and 



v 2 = - , 1 , ,+ 1 



VI + 3e 2( - x+a ^ VI + Se" 2 ^"' 1 ) 
In figures [H and [2J we plot the metric and potential. We see that as a — > oo, the 

metric and potential both tend to 0.5. We expect this since in this limit there is 

negligible interaction between the kinks and there will therefore be geodesic motion 

on Ai with the mass of the system being double the mass of a 6 kink. Indeed the 

asymptotic interaction energy of two <p 6 kinks is 

£ = 2exp(-i?), (19) 

where R is the kink separation. 
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Figure V. The metric g as a function of kink position a. 
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Figure 2: The potential V as a function of kink position a. 



In the next section we consider the dynamical system that is given by the La- 
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grangian (]14p in order to study kink scattering. 



4 The dynamical system 

On varying the Lagrangian (fT4"j) we obtain the field equation 

Ida r, dV . , 

which we can interpret as a particle with position a(t) moving in a potential V(a), 
with variable mass g(a). The initial conditions for kink scattering with a velocity u 
are a(t = 0) = ao and ct(£ = 0) = —it. The total energy of the system is 

E=^g(a )u 2 + V(a ), (21) 

so that the turning point of the motion, a±, is given by V(a\) = E. We can simply 
solve the equation of motion (1201) by quadrature to obtain a(t) implicitly as 




f(a >=/ i/^rW*- < 22) 

which is valid for < t < ti, where ti is the turning time t\ = t(ai). The position 
for t > t\ is determined by the fact that the motion is symmetric about t\, i.e. 
a{t — ti) = a(t 1 — t). 

In order to derive the approximate dynamics, the task is therefore to compute 
the numerical integrals in (122]) . In figure [3] we plot such approximate dynamics for 
ao = 5 and u = 0.3. We also include the dynamics according to a standard Runga- 
Kutta code. We see that the approximate dynamics are qualitatively correct and 
are accurate before the "apparent" point of least approach. The motion according 
to (fi~3]) is seen to take place in a stronger potential than the actual potential. This 
can be compared to Sutcliffe's analysis [10] for the sine-Gordon model where rather 
than using a patching of two kinks like we have done in this chapter, that is 

ip! = 4 arctan e x ~ a + 4 arctan e x+a - 2tt, (23) 
the following patching is used 

tan(^) = e x - a - e~( x+a) . (24) 
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Figure 3: The "real dynamics" (blue curve) for a = 5 and a(t = 0) = —0.3 and the 
"approximate dynamics" (red dots). 

In figure 0J we plot ipi and ^2 with the kinks close together (a = 0.5). We see 
that the kinks, according to ( J24|) . start to deform slightly compared to ( 1231) . This 
has the effect of "softening" the "apparent potential" and is what gives the results 
presented in [10] where both dynamics remain more or less the same for all time. 
The patching for </> 6 kinks presented in this chapter does not, in the same way, take 
account of solitons deforming when they are close to each other. However, the early 
time dynamics according to our approximation show that in this regime it works 
well numerically. Such dynamics have also been calculated for a smaller initial speed 
(u = 0.1) and results obtained are qualitatively the same. 

5 Conclusion 

In this paper, we have investigated the collective coordinate approximation in the 
context of </> 6 kinks with the aim of showing how to test the validity of the method, 
which is much used in the study of solitons. We have found that the Bogomolny 
equations can be solved exactly for this model leading to 1-soliton solutions that 
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Figure 4: Kink fields given by equation (1231) (red) and by equation ( 1241 (green) for 
a = 0.5. 

can then be "patched" together. This field led to a dynamical system for which the 
approximate motion was found via a series of numerical integrals. 

The "real dynamics" were calculated through a numerical code and presented. 
As expected, two kinks moving towards each other at non-relativistic speeds, slow 
down, because of their mutual repulsion, and eventually momentarily stop, then 
move apart. 

With our initial conditions, the approximation we have made leads to dynamics 
which, for time before t « 10, are of a good level of accuracy. Due to the way 
we patch the kinks together, the solitons observe a potential that is stronger than 
the actual potential. We conclude that the collective coordinate approximation 
presented in this paper leads to dynamics that are qualitatively the same as the 
actual dynamics and works well numerically in the early time regime, and that the 
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approximation could be improved if the deformation of solitons as they begin to 
coalesce could be taken into account. 
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